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I review how the phenomenology of localization applies to fermions in lattice gauge theory, present measure- 
ments of the localization length and other quantities, and discuss the consequences for things like the overlap 
kernel. 
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Localization is a phenomenon long studied in 
condensed-matter physics (Ij. Concerned with 
conduction in disordered media, it is based on the 
study of eigenstates of the Schrodinger equation 
in a random potential. Its counterpart in lattice 
gauge theory is the study of the spectrum and 
eigenstates of a Hermitian fermion kernel in the 
fluctuating gauge field. Here we deal with Wilson 
fermions, with kernel Hyy = "/zDyy, in an assort- 
ment of gauge ensembles \2V6\ . We find low- lying 
localized states, extending up to an energy called 
the mobility edge; and extended states at higher 
energy. We characterize the localized states by 
their localization length and their support length 
(related to their inverse participation ratio, or 
IPR). In this talk I will expand a bit on our 
methods for calculating these quantities, which 
are based on calculations of fermionic Green func- 
tions. 

The Wilson fermion kernel is not so popular 
these days for its own sake but rather for its role 
in constructing domain-wall fermions 0] and the 
closely related overlap fermions . The simplest 
overlap kernel is given by 
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1 -75Sgn(i/vy). 
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The range of Dqv can become long if the 
spectrum of Hw contains low-lying (or zero) 
modes Whether this happens, however, de- 
pends on whether the low-lying modes are local- 
ized. Localization is reflected in the Green func- 
tions of the theory. On the one hand, this enables 
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us to use Green functions to calculate localiza- 
tion quantities; on the other hand, the range of 
the overlap kernel, which is just a Green func- 
tion, will be influenced directly by these quanti- 
ties. We will argue here that, though the energy 
threshhold of Hw lie at zero, the range of Doy is 
determined by the position of the mobility edge — 
that is, the low-lying modes are generally harm- 
less. The mobility edge perforce lies at non-zero 
energy, as long as one avoids the Aoki phase |7]. 
Indeed, the descent of the mobility edge to zero 
serves as a useful definition of the onset of the 
Aoki phase. 

1. Localization basics 

The Wilson fermion kernel is 



If (W + amo) -C 



Here 



(2) 
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is the naive Dirac term while 

is the Wilson term. A Hermitian operator can be 
obtained by defining 



Hw — Dwlf5 — h}^ 



(5) 



and it is clear that Hw is a "Hamiltonian" of a 
fermion moving in the random gauge field Ux^- 
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Figure 1. How the spectrum of Hw changes as one adds (top to bottom) one dislocation, many disloca- 
tions, and a random ensemble of dislocations to the gauge field. 



Dislocations in any given configuration U^^ can 
create bound states, and even zero modes 8 , in 
the spectrum of Hyy. Of course we do not study 
single configurations of Uxfi but averages over an 
ensemble. The result is localization. 

For a general overview, let's look at Fig. □ 
which shows how localization can appear as we 
go from free fermions (i.e., fixing Uxfj. = 1) to the 
fully fluctuating gauge field. We assume that the 
free operator Hw has a gap for the parameters 
we choose; above that gap lie plane- wave states. 
Condensed-matter physicists call the top of the 
gap the band edge. If we add a single dislocation 
to the otherwise constant gauge field, the plane 
waves will be replaced by scattering states and 
a bound state might appear in the gap. A large 
number of dislocations will create a large num- 
ber of bound states, and the scattering states will 
show effects of multiple scattering. In finite vol- 
ume, the scattering states as well as the bound 



states form a discrete spectrum. 

In the infinite volume limit, with a fixed den- 
sity of dislocations, both the bound states and 
the scattering states form a continuum. The en- 
ergy that separates the two is called the mobil- 
ity edge. The infinite volume limit automatically 
gives us an average over the shapes and positions 
of the dislocations in any finite subvolume. Al- 
ternatively, we can keep the volume finite and 
average over gauge field configurations ourselves, 
as is done in lattice gauge theory. 

Hw can only have zero eigenvalues in the su- 
percritical region, —8 < mp < 0. As it happens, 
this includes the region of interest for the con- 
struction of domain- wall and overlap fermions. 
The Aoki phase lies in this supercritical region; 
we stay outside the Aoki phase so that the mobil- 
ity edge is above zero. We will present numbers 
below for toq = —1.5. For the gauge couplings 
we choose, the theory lies between the Aoki "fin- 
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gers." These couplings are in fact popular among 
users of domain- wall and overlap fermions. 

2. The spectral density 

Since averaging over the gauge field smears the 
energy eigenvalues, there is never a discrete spec- 
trum of bound states accompanied by a contin- 
uum of extended states. The spectrum in the 
disordered system is described by a continuous 
density of states, 

n 

where A is the energy eigenvalue. We show in 
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Figure 2. Eigenvalue density p(A), in lattice units, 
for f3 = 5.7 (a~^ = 1 GeV, upper curve) and 
(3 = 6.0 (a~^ = 2 GeV, lower curve). Plaquette 
action. 



Fig.[21the density of states for the gauge ensemble 
generated with the Wilson plaquette action. The 
vertical bars indicate the mobility edge for the 
two couplings. Nothing special happens at the 
mobility edge. 

[We find it convenient to calculate p via the 
identity 

TrplA) = -J- lim (ImTrG(A4-imi)), (7) 

V mi— ►O 

where G(z) = {Hiy — z)~^ is the resolvent oi H^. 
The mobility edge is where the averaged localiza- 
tion length li{X) (see below) diverges.] 



Table 1 

Eigenvalue density p(0) at A = for three gauge 
actions, each at two values of the gauge coupling. 



Cutoff a-' 


Wilson 


Iwasaki 


DBW2 


1 GeV 


10-^ 


6 X IQ-^ 


10-^ 


2 GeV 


10-4 


7 X 10-^ 


< 10-^ 



The main point of Fig.|21is to show that there is 
indeed a nonzero density of states at zero energy. 
Table m shows how this density changes with cut- 
off and with the choice of gauge action. Improved 
actions dramatically lower p{0), but in no case is 
it actually zero. 

In a theory with dynamical Wilson fermions, 
these zero modes will be of no consequence be- 
cause they will cause the fermion determinant 
to vanish. In any other theory, however, the 
spectrum of Hyy does not directly determine the 
fermion determinant. Whether the ensemble is 
quenched, or whether it contains the domain- 
wall or overlap fermion determinant, the ensem- 
ble average will furnish Hw with modes near (and 
at) zero energy. Will this destroy the locality 
of Doy? No — because (outside the Aoki phase) 
these modes are localized.^ 

3. Localized modes 

A localized mode can be described by its sup- 
port length Is , which is the size of the region that 
contains most of the mode's density; and by its 
localization length 1^, which is the decay length 
of any long-range tail (see Fig. |3J|. Each can be 
averaged over all eigenmodes with eigenvalue A, 
giving ^^(A) and ls{X). 

The average support length can be used to de- 
termine in what range of A the localized modes 
are dilute. In view of the mode density shown 
in Fig. 12 it is clear that low-lying modes will 
be dilute while the higher modes will be more 
crowded. The average localization length, on the 
other hand, will directly affect the decay rate of 

^This is not to say that these modes won't cause any trou- 
ble at all. Low-lying modes give Hw a bad condition num- 
ber, and make it difficult to calculate the overlap kernel 
via various approximations to Eq. Q- See for instance 
Refs. MM . 
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Figure 3. Characteristic lengths of a locahzcd cigenfunction: the support length Ig, which is the size of 
the region containing most of the amplitude; and the localization length Ig, which is the decay length of 
the tail. 



Green functions of Hw- We turn ^ to these 
Green functions to calculate both li{X) and ls{X). 

We define a "meson" Green function as a pair 
propagator, 

where the complex variable z = X + imi . Intro- 
duction of A will enable us to focus on the contri- 
butions of eigenfunctions with eigenvalue near A; 
the imaginary part mi is a regulator that is put 
in to avoid the singularities on the real axis due 
to these eigenfunctions. 

The spectral representation of F is a double 
sum, 

^<-(^^*"^(^)A;:7-i^)- 

(9) 

In the limit mi ^ it diverges as 1/mi; the 
divergence comes from the terms where n"*" — , 
and thus we are left with a single sum. 




*„(a;)p|^'„(2/) 



The mi — > limit has thus brought out the 
eigenmode densities |5'„(a;)|^. The limit also fo- 
cuses on modes with A„ w A, since the fraction 
within the bracket approaches a delta function as 
mi 0. 

Now we set x = 0, y — (y,t), and sum over 
the spatial coordinate y. This gives the correla- 
tion function F(i) between spatial slices. Taking 
finally the mi ^ limit, we find 

lim miF(t) - 7rp(A) / dlVx{l)e-*/\ (11) 

mi— >0 J 

where 'P\{1) is the distribution of localization 
lengths among the modes with eigenvalue A. The 
tail of T{t) thus gives an average of the tails of 
the mode densities |\E'„(x)p, and the average 
can serve as an estimate of the localization length 
li{X). By definition, ^^(A) diverges at A = Ac, 
the mobility edge. Clearly this happens when the 
modes at A become extended, meaning they have 
spread across the lattice. 

The same tools give us an estimate of the aver- 
age support length ^s(A). If a given (normalized) 
mode density |^„(a;)p takes a value « 1/Zf over 
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a region of linear size 1^, then 



which approaches 



nil 



V\^h{n) (A„-A)2 + m2 
+0(mi), (12) 

npiX) 



as mi 0. Thus we obtain an estimate of 
the average support length. Comparing ls{X) to 
p{X)~^^^, which is the mean distance between lo- 
calized modes, tells us whether the modes at A 
overlap with each other or whether they are di- 
lute. 

4. Back to the overlap kernel 

The dilute, the less dilute, and the extended 
modes all contribute to the overlap kernel -Dqv- 
The contribution of the dilute localized modes, 
those below some energy A, is of the form 



{\Dov{x,y)\)i, 



oc 



dXp{X) 



X exp 



Ap(A) exp 



\x - y\ 
' 2h{X) 

\x -v\ \ 
2h{X) ) ' 



(13) 



where the second line follows from the steep rise 
in p(A) seen in Fig. El Note that these isolated 
modes contribute an exponential tail character- 
ized by their localization length; their low energy 
is harmless. The dense modes are up near the 
mobility edge Ac; if we lump them with the ex- 
tended modes then we can estimate that the fall- 
off of their contribution will be governed by Ac, 



(|i:>ov(a;,y)|)ext «Cexp(-Ac|x-y|), 



(14) 



where C is 0(1). 

Let's put in numbers for the plaquette action 
at (i = 6.0. We measure the mobility edge to be 



at Ac — 0.41 and we fix the demarcation point to 
be A 0.2. The result is 



{\Doy{x,y)\) 



10 exp ( ^ 



+0{l) ■ exp 



2.4 



.(15) 



The second term, governed by the mobility edge, 
wins — both in prefactor and in exponent. We find 
this to be true in all the cases (actions and cou- 
plings) we have studied, including those in Table 
1. The reverse could be true in other cases; then 
the range of -Dov would be determined by the 
representative localization length, li{X). The fact 
that the gap is zero would still be harmless. 

We can interpret Ac as a mass scale, the mass 
of effective excitations that influence the overlap 
kernel. When the cutoff is a^^ = 2 GeV, we find 
that Acfl^^ ~ 800 MeV for all three actions. It 
may be a disappointment to find that the cutoff is 
felt at an energy that is so low, but this is still well 
above Aqcd which is the true dynamical scale. 
One ought to worry, however, at stronger cou- 
phngs. When — 1 GeV, the value of AcC"^ 
comes out to be only 250-320 MeV, not a place 
we would like to see unphysical particles in the 
spectrum. 

5. Conclusions 

The three main lessons of our work: 

1. The mobility edge at Ac > [or le{X), if it is 
greater than A~^] assures a finite range for 
-Dov I even though Hw has no gap in the 
disordered gauge field. 

2. One should demand that XcQ^^ ^ Aqcd- 

3. AcC^^ is fairly insensitive to the gauge ac- 
tion (at fixed cutoff), even as p{0) varies 
widely. 

The third point emerges from numerical results 
|2| that are not presented here; it implies, ac- 
cording to our arguments, that the range of the 
overlap kernel will be equally insensitive to the 
gauge action. This has already been reported 
in Ref. [HI; moreover, the range of Dov given 
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there for w 2 GeV agrees well with our es- 
timate, shown in Eq. (|15|l above. It is curious 
that the range given for the overlap kernel when 
fti 1 GeV in Ref. m] is much shorter (by a 
factor of 4) than our estimate would indicate. The 
reason for this is an interesting open question. 

I have not discussed the implications of local- 
ization for domain-wall fermions, where it would 
govern both the range of the equivalent overlap- 
like kernel and the residual mass due to finite- 
ness of the fifth dimension in practical calcula- 
tions |2i;ij . Numerical results from QCDOC were 
reported by Peter Boyle at Lattice 2005 (12j . 
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